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ABSTRACT 

We explore the effects of active galactic nuclei (AGN) and star formation activity on the infrared 
(0.3 - 1000 fim) spectral energy distributions of luminous infrared galaxies from z = 0.5 to 4.0. We 
have compiled a large sample of 151 galaxies selected at 24 fim (5 2 4 ^ 100 ^iJy) in the GOODS-N 
and ECDFS fields for which we have deep Spitzer IRS spectroscopy, allowing us to decompose the 
mid-IR spectrum into contributions from star formation and AGN activity. A significant portion 
(<~ 25%) of our sample is dominated by an AGN (> 50% of mid-IR luminosity) in the mid-IR. 
Based on the mid-IR classification, we divide our full sample into four sub-samples: z ~ 1 star- 
forming (SF) sources; z ~ 2 SF sources; AGN with clear 9.7 /im silicate absorption; and AGN with 
featureless mid-IR spectra. From our large spectroscopic sample and wealth of multi-wavelength data, 
including deep Herschel imaging at 100, 160, 250, 350, and 500 fim, we use 95 galaxies with complete 
spectral coverage to create a composite spectral energy distribution (SED) for each sub-sample. We 
then fit a two-temperature component modified blackbody to the SEDs. We find that the IR SEDs 
have similar cold dust temperatures, regardless of the mid-IR power source, but display a marked 
difference in the warmer dust temperatures. We calculate the average effective temperature of the 
dust in each sub-sample and find a significant (~ 20 K) difference between the SF and AGN systems. 
We compare our composite SEDs to local templates and find that local templates do not accurately 
reproduce the mid-IR features and dust temperatures of our high redshift systems. High redshift IR 
luminous galaxies contain significantly more cool dust than their local counterparts. We find that a 
full suite of photometry spanning the IR peak is necessary to accurately account for the dominant 
dust temperature components in high redshift IR luminous galaxies. 

Subject headings: Galaxies: evolution - Galaxies: active - Galaxies: star formation - Infrared: galaxies 
- ISM: dust 
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1. INTRODUCTION 

Understanding the link between star formation and 
active galactic nucleus (AGN) activity in the high red- 
shift Universe is key to determining a complete picture 
of galaxy evolution. The history of the star formation 
rate density i n the Universe show s a peak in the range 
1-3 (e.g., Bouwens et al. 2009 1, during which time at 



least 50% of the stars in th e local Universe are ex pected 
to have been formed (e.g., Dickinson et al. 2003). Fur- 



thermore, the bulk of the star tormation during this peak 
period is occurring in luminous infrared galaxies (LIRGs, 
L m = 10 n -10 12 L Q ) and ultra lu minous infrared galax- 
ies (ULIRGs, L T r > 10 12 L Q ) (e.g., |Magnelli et al.(|2011 
Murphy et al.| 2011 ) . During the same epoch, the black 
holes within tr ie centers of massive galaxies are building 



up their mass fWall et al. |2005| [Kelly et al.|2010| , making 



this a key period for the growth of galaxies 

It is important to understand how each mechanism 
(star formation and AGN activity) manifests itself in 
observations of the interstellar medium, particularly the 
dust. In the far- infrared, dust grains emit thermal emis- 
sion characterized by a blackbody spectrum with an ad- 
ditio nal v$ term to ac count for the emissivity of the 
dust ( |Hildebrand| [l983 ) . Each galaxy will have a distri- 
bution of dust temperatures depending on the size and 
distribution of the dust relative to the heating sources; 
the integrated IR spectral energy distribution (SED) will 
be a combination of modified blackbodies from all dust 
components. Most studies of high redshift galaxies fit a 
single modified blackbody template, assuming an aver- 
age dust temperature, but there is evidence from local 
star-forming galaxies and ULIRGs that at least two dust 



SED (Dunne & Eales 200llparrah et al. 


2003 


Willmer 


et al. 


2009( |Galametz et al.pDTl|). Phys 


ically, this can 



be found in a galaxy: cool diffuse (cirrus) dust in the 
interstellar medium; warmer dust in active star-forming 
regions; and even warmer dust found in an AGN torus 



( Rowan- Robinson & Crawford 1989). 

'lb study the properties of AGN and star-forming (SF) 
galaxies at high redshift, it is necessary to first identify 
systems likely harboring an AGN. Many previous stud- 
ies have used X-ray detections to determine whether a 



galaxy is hosting an AGN (e.g., Alexander et al.||2003 [) 
However, as the majority of star tormation at z ~ 1-3 
occurs behind dust, many AGN will also be obs cured by 
dust and gas an d undetectable in the X-ray (e.g., Alexan- 



der et al. 2008). Therefore, we must rely on alternative 
methods which are insensitive to dust obscuration, to 
identify the presence and strength of an AGN. 

The mid-infrared spectrum, consisting of dust- 
reprocessed light, is rich in information about the under- 
lying power sources of the galaxy. The mid-IR spectrum 
contains three major features that can be used as a di- 
agnostic of the mechanism driving the intense luminos- 
ity of the galaxy. First, the mid-IR spectrum exhibits 
polycyclic aromatic hydrocarbon (PAH) emission lines. 
PAHs are excited by UV and/or optical photons, most 
often produced by young stars; as such, PAHs are good 
tracers of the amount of star formatio n in a galaxy (e.g., 
Peeters et al.|2004||Brandl et al.|20"06| ). Second, the mid- 
1R spectrum may exhibit additional continuum emission 



comin g from hot dust around the AGN. [Mullancy et al.| 
(20111 find that, on average, such emission can be rep- 
resented as a steep power law {v~ 2 ) below 20 /im which 
flattens at longer wavelengths, falling off steeply beyond 
40 /im. The final main features observed in the mid-IR 
spectrum of galaxies are due to silicate absorption, most 
clearly seen at 9.7 jum. This absorption feature can be 
due to either the AGN or the star formation, since it is 
produced by a dust screen surrounding a hot emission 
region, though there is some observational evidence that 
it is m ore prevalent in AG N with an optically thick dusty 
torus flSpoon et aL]|2007| ). 

High redshift (z = 1 — 3) (U)LIRGs are an attractive 
option for studying both SF and AGN activity in this 
peak epoch of stellar growth. Not only are they undergo- 
ing intense periods of star formation (SFR > 10 - 100 M Q 
yr -1 ), but many show signs of concurrent A GN growth 



(e.g., Lutz et al. 2008 Coppin et al. 20101. With the 



advent of the space-based IR observatories Spitzer Space 



Teles cope and Herschel Space Observatory (Pilbratt et al 



20101, we are now in a position to obtain spectroscopy 
and photometry for these systems from near- to far-IR 
wavelengths, which will allow us to decompose the SEDs 
into the relative contributions from AGN and SF activity. 
Quantifying each component in high redshift (U)LIRGs 
will enable us to investigate evolutionary trends of SF 
galaxies and AGN with redshift and determine how the 
two components might be affecting the observed dust 
properties 

We have assembled a sample of 151 high redshift 
(U)LIRGs with deep Spitzer mid-IR spectroscopy. We 
use the mid-IR spectrum not only to determine the domi- 
nant mid-IR power source (SF or AGN activity), but also 
to quantify the relative contribution of each to the mid- 
IR luminosity of every galaxy. We combine our spec- 
troscopy with multi-wavelength data, specifically new 
Herschel imaging at 100, 160, 250, 350, and 500 /im, 
to study the full IR SED of high redshift LIRGs and 
ULIRGs. We aim to study how mid-IR classification is 
related to far-IR properties such as dust temperature and 
IR luminosity. 

High redshift (U)LIRGs have long been taken to simply 
be analogs to local systems of the same luminosity. Due 
to the limited data previously available for high redshift 
(U)LIRGs, astronomers often resort to blindly applying 
local galaxy templates to high redshift systems to esti- 
mate parameters such as dust temperature and total IR 
luminosity. Recent studies have found that the SEDs of 



at z > 1.5 (e.g., 


F 


'ope et al.||2006| iPapovich et al.||2007l 
Dannerbauer et al.||2010| |Elbaz et al.| 


Desai et al.||200/ 


2010J Noruon et a 


. 2010). Given the considerable size of 



of multi- wavelength data, we are in a unique position to 
create full IR SEDs for high redshift (U)LIRGs, test how 
accurately they reflect what we see in the local Universe, 
and determine where and when variations occur. 

The paper is laid out as follows. In Section 2, we 
describe details of our sample and the mid-IR spectral 
decomposition of individual sources. In Section 3, we 
describe how we created composite IR SEDs for four 
sub-samples. In Section 4, we discuss specific features 
of our new composite SEDs, and we compare these com- 
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posites to local templates. Finally, in Section 5 we dis- 
cuss our results and summarize in Section 6. Through- 
out this paper, we assume a standard cosmology with 
H = TOkms- 1 Mpc^ 1 , fi M = 0.3, and Q A = 0.7. 

2. DATA 
2.1. Sample Selection 

Our sample consists of 151 high redshift galaxies from 
the Great Observatories Origins Deep Survey North 
(GOODS-N) and Extended Chandra Deep Field Survey 
(ECDFS) fields. We include all sources in these fields 
that were observed with the Spitzer I RS (for details on 
the instrument, see |Houck et al.1|2004| . While this sam- 
ple contains a diverse range of sources, depending on the 
goals of each individual observing program, the overly- 
ing selection criterion is that each source is bright enough 
at 24 /im to be observable with the IRS in a reasonable 
amount of time (< 10 hours). As a result, 93% of the 
sources have S24 > 100 jttJy, 80% have S24 > 200 /iJy and 
60% have S24 > 300 /iJy. More details on this database 
of IRS sources in GOODS-N and ECDFS can be found 
in a separate paper (Pope et al. in preparation). 

It is important to determine how representative the 
IRS sample is of the full MIPS 24 /im population. Within 
GOODS-N, the IRS sources make up ~ 20% of all 24 /im 
sources above 200 /zJy and ~ 30% above 300 /iJy. If we 
limit ourselves to only 24 /im sources with spectroscopic 
redshifts greater than 0.5, then the IRS sources make 
up ~ 40% above 200 /iJy and ~ 60% above 300 ,uJy. In 
Figure [T] we show that our IRS sample spans the same 
distribution in Herschel+ Spitzer colors as the full MIPS 
population with S24 > 100 /iJy indicating that our IRS 
sample is representative within these parameters. In Sec- 
tion[3j we use 95 of our sources to create composite SEDs, 
and we show the distribution of these template sources 
in Fig. JT] as well. We run a Kolmogorov-Smirnov (KS) 
comparison test of the parent MIPS population and the 
two subpopulations (our full IRS sample and limited tem- 
plate sources). The significance level of the KS test, for 
each of the subpopulations and the MIPS population, is 
93%; therefore, we cannot reject the null hypothesis that 
the populations are drawn from the same parent distri- 
bution. 

As an additional test, we can determine if the frac- 
tion of AGN-dominated sources (defined using the mid- 
IR spectrum, see ^2.3) within our IRS sample is pro- 
portionate to that of all 24 /im sources. In Kirkpatrick 
et al. (2012b, ApJ submitted), we develop new color- 
color diagnostics to separate AGN and SF sources at 
high redshift calibrated using this IRS sample. Using 
our new color cuts, we determine which 24 /zm sources 
from the full MIPS samples are AGN dominated in the 
mid-IR. The AGN dominated sources in our IRS sample 
have a flux density distribution consistent with the mid- 
IR color-selected AGN in the full MIPS sample above 
5*24 > 100 jiiJy. This gives us confidence that the IRS 
sample is providing fair and unbiased sampling of the 
full 24 /im population above these flux limits. We report 
the Spitzer MIPS 24 and 70 /jm photometry and Her- 
schel PACS and SPIRE photometry for our full sample 
in Appendix A. 

Furthermore, it is important to determine the overlap 
between 24 /zm and 100 /im selected galaxy samples to de- 



termine how representative our 24 /im selected sample is 
of the greater population of IR luminous galaxies at high 
redshift. When looking at a blind catalog of sources se- 
lected at 100 /im in GOODS-N with S wo > 1.1 mJy (the 
detection limit for this survey, see Elbaz et al. 2011), we 
find tha t 97% of PACS 100 /im are detected at 24 /im 



(see also |Magdis et7^|2lHlj ) , and 67% of PACS 100 /im 
sources have S24 > 101) /iJy. Moreover, 70% of our IRS 
sources have a detection at 100 /mi. Therefore, while 
we are selecting sources at 24 /im in this study with 
S*24 > 100 /iJy, we are getting a representative sampling 
of the bulk of the PACS 100 /im selected sources. Given 
the large beam sizes at SPIRE wavelengths, which make 
robust counterpart identification challenging, we are un- 
able to perform a similar simple test of the overlap of 
24 /im sources with those selected in SPIRE. However, we 
do note that the bulk of submm sources selected at even 
longer wavelengths, 850 /im, are detected at 24 /im and 
most of these det ections are brighter than S24 = 100 /iJy 
( Pope et al. 2006 ) . It is important to keep in mind that in 
this paper, we focus on sources that have mid-IR spec- 
tra and PACS and/or SPIRE photometry. Therefore, 
our sample is representative of sources that are detected 
both in the mid-IR and far-IR and may not cover the 
parameter space of sources fainter than our flux limits or 
sources detected in either the mid-IR or the far-IR. 

The low resolution (R = A/AA - 100) Spitzer IRS 
spectra were reduc ed following the method outlined in 
Pope et al. (2008a). Specifically, since many of these are 
long integrations, we take care to remove latent build- 
up on the arrays over time, and we create a 'supersky' 
from all the off-nod observations to remove the sky back- 
ground. One dimensional spectra are extracted using the 
Spitzer IRS Custom Extraction (SPICE) in optimal ex- 
traction mode. For each target a sky spectrum is also 
extracted to represent the uncertainty in the final target 
spectrum. 

2.2. Multi-wavelength Data 

The GOODS fields have been extensively surveyed and 
are rich in deep multi-wavele ngth data, including: Chan- 
dra 2 Ms X- ray observa tions ( [Alexander et al.|2003 Luo 
eFld71[2008j |Xue et al.|[20TT| ); 3.6, 4.5, 5.8, and WXTpm 
imaging from the Infrared Array Camera (IRAC) on Sp- 
tizer (Dickinson et al. in preparation); IRS peak-up ob- 
servations at 16 /im ( |Teplitz et al.|20il|); and MIPS imag- 
ing at 24 and 70 /im (|Magnelli et al.||201l[ ). Recently, 
GOODS-N and GOODS-S have been surveyed with the 
GOO DS- H erschel Open Time Key Program ( P.I. David 
Elbaz ,|Elbaz et al.|2011[ ) using both the PA CS ( |Poglitsch 
et al.|2010p and SPIRE ( jGriffin et al.|2010[ ) instruments 
providing deep photometry at five far-IR and submm 
wavelengths: 100, 160, 250, 350, and 500 /im. For the 
Herschel imaging, flux densities and the associated un- 
certainties were obtained by point source fitting using 
24 /im prior positions, allowing us to probe deeper lim- 
its in the Herschel images. In addition, the GOODS- 
Herschel catalog is only comprised of sources with a clean 
detection, based on the 24 /im prior position having no 
bri ght neighbors in a given passband (for further details, 
jElbaz et al.|[20TT| ). 



We combine this space-based imaging with ground- 
based imagin g in the near-IR ( J and K bands) from 
VLT/ISAAC (Retzlaffetal.|2010[ ) and CFHT / WIRC AM 
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Fig. 1. — The distribution of the MIPS 24 /im GOODS sample with S24 > 100 ^Jy and S250 > 3<r and the IRS sample presented in this 
work in a colorspace combining far-, mid-, and near-IR fluxes (left) and in the color S250/S24, (right). In §3] we use 95 sources from our 
IRS sample to create composite templates. We indicate the distribution of these 95 sources in red (see the online journal for a color version 
of this figure). The IRS sample follows the a similar distribution of colors as the full 24 /im population with 524 > lOO^tJy. 



TABLE 1 

Basic properties of our four sub-samples. 



Sub-sample 


No. of 
Sources 51 


Median 
Redshift 


Median S s h 

(MJy) 


Median S 2 4 b 

(MJy) 


Median 5ioo b 
(mjy) 


z ~ 1 SF galaxies 
z ~ 2 SF galaxies 
Silicate AGN 
Featureless AGN 


39 
30 
17 
9 


1.0 [0.8, 1.0] 
1.9 [1.8, 2.1] 
1.9 [1.6, 2.0] 
1.2 [0.6, 1.6] 


42 [23, 57] 
17 [13, 38] 
64 [21, 140] 
288 [240, 311] 


370 [260, 570] 
270 [220, 370] 
470 [250, 860] 
1520 [1240, 2300] 


7.9 [4.8, 14.7] 

3.2 [1.6, 5.0] 

5.3 [3.0, 10.2] 
9.5 [4.0, 11.7] 



We list the lower and upper quartile values in parenthesis next to each median value. 
a We list the number of sources in each sub-sample that are used to create the composite SEDs as 
well as the number of sources we have rejected in each subsample in Section [3] 
b Observed frame flux densities. 



Lin et 
wavelength's 



al. 



wc use 



2012] . At the longest 
available data from 



(|Wang et al.||2010 

(subjmm 

LABOCA on APEX fWeiB et al.|2009| ) and the combined 
AzTE C+MAMBO mm map of GOODS-N flPenner et alT 
2011]. 



2.3. Mid-IR Spectral Decomposition 

We perform spectral decomposition of the mid-IR spec- 
trum (5-12 fj,m rest frame) for each source in order to 
disentangle the AGN and SF com ponents. We follow 



the technique outlined in detail in Pope et al. 

We fit the individua. 



(2008a) 



spectra 



which we summarize here 
with a model comprised of three components: (1) the 
star formation component is represented by ei ther the 
local starburst composite of 
ply th e mid-IR spectrum o 



Brandl et al. ( 2006 ) or sim- 
tfi e pro totypical starburst 



M 82 ( |Forster Schreiber et aT]|2003[ ) - with the SNR, 
wavelength coverage, and spectral resolution of our high 
redshift spectra, both give equally good fits to the SF 
component of our galaxies; (2) the AGN component is 
determined by fitting a pure power-law with the slope 
and normalizat ion as free paramete rs; (3) an extinction 



curve from the Draine et al. (2003) dust models is ap 



plied to the AGJN component. The extinction curve is 
not monotonic in wavelength and contains silicate ab- 



sorption features, the most notable for our wavelength 
range being at 9.7 /im. The local starburst composite 
and the M 82 template already contain some intrinsic ex- 
tinction. We tested applying additional extinction to the 
SF component beyond that inherent in the templates and 
found this to be negligible for all sources. We fit all three 
components simultaneously and integrate under the star- 
burst spectrum and power-law continuum to determine 
the fraction of the mid-IR luminosity (~ 5-12 /im de- 
pending on the redshift of the source) from SF and AGN 
activity, respectively. For each source, we quantify the 
strength of the AGN in terms of the percentage of the 
total mid-IR luminosity coming from the power-law con- 
tinuum component. Based on this mid-IR spectral de- 
composition, we find that 38 (25%) out of our sample of 
151 galaxies are dominated (> 50%) in the mid-IR by 
an AGN. Figure [2] shows examples of the best fit models 
(red dashed line) for an SF-dominated galaxy and two 
AGN-dominated galaxies, with and without prominent 
9.7 /xm extinction. The extincted power law component 
is shown by the blue dashed line, and the PAH template 
is the green dashed line. 

To more thoroughly compare the mid-IR spectral prop- 
erties and far-IR SEDs within our sample, we divide our 
galaxies into four sub-samples based on the results of 
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Fig. 2.— Mid-IR spectra of GNJRS26 (left), GNJRS12 (middle), and GNJRS30 (right). The red dashed curves show the best fit SED 
which is made up of an extincted power law component (blue dashed line) and a starburst template (green dashed line). GNTRS26 has 
strong PAH features indicative of star formation activity, while GN_IRS12 and GN_IRS30 are dominated by the power-law component, 
indicating the presence of an AGN. GNJRS30 shows a strong 9.7/im Si absorption feature on top of the continuum component. The 
relative contributions of each component to the mid-IR luminosity determines whether a galaxy is dominated by AGN or star formation 
activity in the mid-IR. 



the mid-IR spectral decomposition. First, each galaxy 
is classified as either SF- or AGN-dominated, based on 
having < 50% or > 50% AGN contribution to the mid-IR 
luminosity, respectively. We further divide the SF galax- 
ies into two bins: z ~ 1 [z < 1.5); and z ~ 2 (z > 1.5). 
For AGN, the mid-IR spectral features have been pre- 
dicted to reveal the shape of the torus surrounding the 
AGN. A clumpy torus produces a power-law spectrum 
and possibly weak silicate absorption, while the presence 
of stro ng silicate absorption suggests a thick obs curing 



torus (Levenson et al. 



2007 Sirocky et al. 



2008). We 



therefore classify the AGJN according to the shape of their 
mid-IR spectrum; those with measurable 9.7 /im silicate 
absorption (hereafter referred to as silicate AGN), and 
those without (hereafter referred to as featureless AGN) , 
which we have classified by eye. We have a much smaller 
number of AGN sources, so separating further accord- 
ing to redshift would not produce a meaningful sample 
with which to determine the average properties. We are 
unable to classify four AGN sources as they lack spec- 
tral coverage in the relevant range (9-10 /im), so we are 
incapable of determining whether they exhibit silicate 
absorption. We refer to these as unclassifiable AGN in 
the relevant figures. Our four sub-samples are listed in 
Table [I] along with their median redshifts and 8, 24, and 
100 /im flux densities. 

While the majority of our sources that are classified 
as SF-dominated, based on the mid-IR spectra, have a 
negligible (< 20%) contribution from an AGN, the AGN- 
dominated sources exhibit varying degrees of concurrent 
SF activity. Figure [3] shows the distribution of mid-IR 
AGN fraction for the silicate (top panel) and featureless 
AGN (bottom panel). The featureless AGN (lacking sil- 
icate absorption) have a very strong AGN continuum, 
accounting for 80-100% of the mid-IR emission, whereas 
the silicate AGN have a more uniform distribution of 
AGN fraction, with some silicate AGN also having weak 
PAH features. The broad range of AGN fraction for the 
silicate sources is partly a reflection of the difficulty in 
disentangling PAH features from the 9.7 /im silicate ab- 
sorption feature in low SNR spectra. 



10 



o 
co 



o 

CO 



Silicate AGN 




60 70 80 90 
Strength of AGN (%) 



100 



Fig. 3. — Sources whose mid-IR spectrum is dominated by an 
AGN continuum are further divided according to whether they 
exhibit significant 9.7 /im silicate absorption. The majority of the 
sources completely dominated by AGN emission (> 95%) exhibit 
a featureless power-law spectrum (bottom panel) , whereas sources 
with a silicate absorption feature are more uniformly distributed 
according to the strength of the AGN in the mid-IR (top panel). 

2.4. Spectroscopic Redshifts 

We determine redshifts for the majority of our sample 
by fitting the po sitions of the main PAH features (see 
Pope et al. 2008a I . In the case of featureless spectra (only 
9/151 so urces), we adopt available optical spectroscopic 
redshifts ( |Szokoly et al.|2004}|Barger et al.|2008||Popesso| 



G 



et al.|[2009 Stern et al. in preparation), with the excep- 
tion ot one source for which we only have a photometric 
redshift. More details on the optical and IRS-based red- 
shifts of individual sources in our sample are given in 
Pope et al. (in preparation). We compare our redshifts 
derived from the IRS spectrum with the optical spectro- 
scopic redshifts using the parameter Az/(1 + Zirs)- We 
find the mean value to be 0.008 with a standard deviation 
of 0.07. Therefore, we conclude that our IRS redshifts 
are good to ~ 0.07. 

The redshift distribution of our sources is illustrated in 
Figure [4] The redshift distribution exhibits a bimodal- 
ity with two peaks at redshifts of around 1 and 2. The 
bimodality of the reds hift distribution i s a result of the 

At 



2008) 



2, the 



24 fim selection (e.g., Desai et al. 
24 fim observed frame flux density corresponds to the 
8 /im flux density in the source's rest frame. Galaxies 
with intense star formation will have an enhanced 8 fim. 
flux density due to prominent PAH complexes at 7.7 
and 8.6 fim. Similarly at z ~ 1, observed frame wave- 
length 24 /im corresponds to a rest frame wavelength of 
12 fim. The strong PAH emission features at 11.3 /jm 
and 12.7 /zm, present in SF galaxies, will boost the ob- 
served 24 /jm flux density. Fig. [4] shows that our sources 
at the highest redshifts are all dominated by an AGN in 
the mid-IR. While the featureless AGN are distributed 
roughly equally across our redshift range, there is a slight 
dearth of silicate AGN at z ~ 1.5 due to the 9.7 fim ab- 
sorption feature falling into the 24 ^m band. 

In the sections that follow, we use the median redshift 
for each subsample (listed in Table [T]) to compare general 
properties of the subsamples and investigate evolution- 
ary connections between the populations. For our SF 
galaxies, Fig. [4] illustrates that the median redshifts of 
1.0 and 1.9, respectively, are accurate for each subsam- 
ple. The median redshift for the silicate AGN is 1.9, and 
again, the distribution is peaked around this value. How- 
ever, the featureless AGN have a fairly uniform redshift 
distribution with the mean and median redshifts both 
giving z — 1.2. It is important to keep in mind that as 
we did not separate the AGN based on redshift, we are 
not claiming that our subsamples are representative of 
all AGN at each of the median redshifts. 

2.5. Stellar Masses 

We use optical and near-IR photometry to estimate 
the galaxy stellar masses in our sample, and compute 
the median stellar mass for each sub-sample (see Table 
[3]). The stellar masses for the current sample are part of 
a larger catalog of stellar masses, photometric redshifts, 
and multi-wavelength data for galaxies in the GOODS 
fields, a detailed description of which will be given in a 
separate paper (Pannella et al., in preparation). In the 
following, and for the sake of completeness, we will briefly 
describe the sources of the photometry used to calculate 
the stellar masses and the procedures used. 

We have built a PSF-matched multi-wavelength cat- 
alog with 10 passbands from U to 4.5 /im. The opti- 
cal, near-IR, and IRAC data used to calc ulate the stel- 
lar masses in GOODS-N are presented in |Capak et al.| 



(2004); 



Wang et al.| pOlOj ); |Lin et alT| j2012p ; and Dick 
mson et al. (in preparation), respectively. A i\ " s -band 
selected ca talog h as been built using SExtractor (jBertin 
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Fig. 4. — Redshift distribution of sources. The distribution ap- 
pears to be bimodal with peaks at z ~ 1 and z ~ 2 reflecting 
the 24 fim selection criterion. At redshifts 1 and 2, the observed 
frame 24 fim band corresponds to the rest-frame 12 fim and 8 fim 
band, respectively. In SF galaxies, those wavelengths are domi- 
nated by prominent PAH features if they lie at these particular 
redshifts, causing them to appear more luminous at 24 fim. The 
sources have been separated according to four sub-samples based 
on the mechanism driving the mid-IR luminosity. Our three high- 
est redshift sources are all AGN, though due to lack of coverage at 
9.7 fim, we are unable to classify them. 



the procedure describe d in Pannella et al. (2009b) and 



Strazzullo et al. (2010). For the ECDFS sources we have 
used photoinetx y from both the GOODS-M USIC ([San 



I tini et al.l[2009| and the GOODS-MUSYC ( |Cardamone 
|et al . 2010) surveys to estimate stellar masses. 

Stellar masses are estimated for bo th fields using the 
SEP fitting code described in detail in Drory et al. (2004, 
20091. We parametrize the possible star formation his- 



tories by a two-component model, consisting of a main, 
smooth component described by an exponentially declin- 
ing star formation rate (tp(t) cx exp(t/r)), linearly com- 
bined with a secondary burst of star formation. The 
main component time-scale r varies in G [0.1, 20] Gyr, 
and its metallicity is fixed to solar. The age of the main 
component, t, is allowed to vary between 0.01 Gyr and 
the age of the Universe at each object's redshift. The 
secondary burst of star formation, which cannot con- 
tain more than 10% of the galaxy's total stellar mass, 
is modelled as a 100 Myr old constant star formation 
rate e pisode of solar metallicity. We adopt a |Salpeter"| 
( 1955 1 initial mass function for both components, with 



& Arnouts 1996) in dual image mode, closely following 



lower and upper ma ss cutoffs of 0.1 and 1 00 M Q , respec- 
tively. Adopting the Calzetti et al. ( 2000 ) extinction law, 
both the main component and the burst are allowed to 
exhibit a variable amount of attenuation by dust with 
Ay € [0, 1.5] and [0, 2] for the main component and 
the burst, respectively. We confirmed our masses by also 
fitting without the secondary burst component and find 
consistent results. 

A significant fraction (~ 66%) of our sources that are 
AGN-dominated in the mid-IR do not display a promi- 
nent stellar bump. For these objects, we do not expect 
the photometry to diverge from the best-fit results until 
after 1.5 fim, at which point the stellar bump becomes 
diluted. Based on our fitting method, we find that this 
should produce a bad fit rather than a biased one. We 
test our fitting method for power law sources on an indi- 
vidual basis by removing the IRAC photometry to ensure 
we were only fitting wavelengths bluer than the 1.6 fim 
bump. The results are consistent with the results pro- 
duced when the IRAC data is included. A recent study 
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by Hainline et al. (2012) found that for IRAC power-law 
sources, the stellar mass can be a factor of 1.4 higher 
when an AGN component is not accounted for in the 
SED modeling. In this study, we use the stellar masses to 
calculate specific star formation rates to as sess w hether 
our sources lie on the main sequence (see f 5.2.1). The 
possible bias of 1.4 introduced by not including an AGN 
component in the stellar fitting will not greatly affect the 
discussion in that section given the spread in the main 
sequence. 

In Section [5] we com pare with stellar masses presented 
(20111. Our stellar masses were calcu- 



Elbaz et a. 



lated using the stellar pop ulatio n synthesis models f rom 



Bruzual fc Charlotl (|2003| , w hile |Elbaz et al | J2011[) use 



templates from pEGaSE.2 (Fioc & Rocca-Volmerange 
19991. As a result, our m asses are a factor ot ~ 2 lower 
(see also |Bell et al.|[20TT| ) . 



3. COMPOSITE SPECTRAL ENERGY DISTRIBUTIONS 

As described in detail above, we divide our full sample 
into the four sub-samples listed in Table [I] For each sub- 
sample, we create a composite SED from 0.3 to 600 /im 
rest-frame by combining data from: ground-based near- 
IR; Spitzer IRAC, IRS, and MIPS (24 /jm, 70 /jm); Her- 
schel PACS and SPIRE; and ground-based 870 /im and 
1.15 mm. 

For sources lacking a detection at the Herschel wave- 
lengths, we extract a measurement of the flux density 
and associated uncertainty for each source directly from 
our images using the prior position at 24 /im. This is not 
just an upper limit but rather an actual best estimate of 
the flux of the source even if it not formally detected. We 
take measurements from the image to avoid biasing our 
SEDs high in the submm by including only detections. 
We do not take a measurement when a source looks too 
blended on the image; we have rejected 26 sources from 
our composite SEDs due to blending in the Herschel im- 
ages. 

Since spectroscopic redshifts are known for all sources 
in our sample, we begin by shifting all spectra and pho- 
tometry to the rest frame. We then examine the full 
IR SED of each galaxy; sources which lack a sufficient 
amount of data in the far-IR/submm (< 3 measurements 
beyond a rest wavelength of 25 /jm) are not used in cre- 
ating the composite SEDs to avoid biasing our results. 
By requiring sources to have measurements (not neces- 
sarily detections) at longer wavelengths, we have only re- 
jected a limited number of sources (26, or 17%) which are 
severely blended with other sources in the large SPIRE 
beams. Our rejected sources have the same redshift dis- 
tribution as our full sub-samples, ensuring that we have 
not introduced a redshift bias. 

When combining the spectra and photometry for the 
composite SEDs, we normalize the data for each source 
to the mid-IR spectrum over the wavelength range 6.4 
to 7.5 /im (with the exception of the featureless AGN), 
chosen because it is free of prominent PAH features, and 
because the majority of our sources have spectral cover- 
age in this range. Any sources with exceptionally noisy 
data over this range (such that we are unable to reli- 
ably decompose the AGN and SF components) or with- 
out continuous data from 6.4 - 7.5 /im are excluded from 
the composite SEDs. The featureless AGN have mid- 
IR spectra that are free of PAH and silicate features, so 



we choose to normalize over a longer wavelength range, 
6.4 - 11.4 /im. 

For clarity, we now state exactly how many sources are 
rejected from the creation of the composites in each sub- 
sample. In the z ~ 1 SF galaxies subsample, we begin 
with 69 sources; we reject 13 due to blending in the far- 
IR and 17 that lack coverage in the mid-IR normalization 
range; this leaves 39 sources. For the z ~ 2 SF galax- 
ies, we begin with 44 sources; we reject 7 sources due to 
far-IR blending, 5 sources that lack coverage in the mid- 
IR normalization range, and 2 sources with exceptionally 
noisy mid-IR spectra; this leaves 30 sources. For the sil- 
icate AGN, we begin with 22 sources; we reject 4 due to 
blending in the far-IR and 1 because it lacks spectral cov- 
erage in the mid-IR normalization range, leaving us with 
17 sources. Finally, in the case of the featureless AGN, 
we initially have 12 sources, from which we reject two 
that are blended in the far-IR and one because it lacks 
coverage in the mid-IR normalization range; we therefore 
have 9 sources which we use to create the composite. The 
number of sources comprising each composite are listed 
in Table [T] The remaining 95 sources are representative 
of our full IRS sample: 82% have S 2 4 > 200 /iJy, and 
60% have S24 > 300 //Jy. For all four subsamples, the 
redshift distributions are the same for the full subsample 
and the limited subsample used to create the composite. 
The specific sources we have rejected in each subsample 
are indicated in Appendix A. 

We determine the median L u for each source over this 
wavelength range. Within each sub-sample, we also de- 
termine the median L u for all sources, and we normalize 
the rest-frame data of each source to the median L v of 
that sub-sample in order to preserve the intrinsic average 
luminosity of each sub-sample. 

In creating the composite SEDs, wc treat the 
near/mid-IR data separately from the far-IR/submm 
data. This is because in the far-IR the individual photo- 
metric uncertainties dominate the noise, whereas in the 
near/mid-IR the scatter between different sources in each 
composite is a larger source of error. After normalization, 
we average the near-IR and mid-IR data by determin- 
ing the median L v and wavelength in bin sizes > 0.1 /im 
(larger than the spectral resolution of the IRS spectra), 
chosen so that each bin is well-populated (> 10 data 
points). We bootstrap the sample to estimate the errors 
on our composite in the mid-IR and near-IR, since in this 
regime the uncertainty in the composites is dominated 
by the scatter between different sources. For each sub- 
sample, we randomly draw sources with replacement and 
recalculate the normalized median SED 10,000 times. 
Because we normalize in the range 6.4 — 7.5 /im, the re- 
sulting SEDs exhibit little scatter in and around these 
wavelengths. 

The far-IR data has two sources of error: the instru- 
ment noise and the confusion noise. We account for both 
in our fitting by adding the instru ment noise in quadra - 
ture with the confusion noise (from Nguyen et al.| ( |2010[ ) ) 
for each flux measurement. We fit all of the normalized 
far-IR (> 18 /im) data in each sub-sample with a two 
temperature component modified blackbody function of 
the form 

M v =aix/x B v (T waIm ) + a 2 x v p x 5„(T cold ) (1) 
where B u is the Planck function. We keep the emis- 
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sivity fixed at f3 — 1.5. Since the photometric uncer- 
tainty is large in the far-IR, we use a Monte Carlo sim- 
ulation to determine the parameters and errors of the 
blackbody model. We randomly sample each data point 
within its error and refit the model 10,000 times. We 
fit for the normalization factors a\ and together with 
T warm and T" C oid , the temperatures of the warm and cold 
dust components, simultaneously using a \ 2 minimiza- 
tion technique. The final parameters for each compos- 
ite are the mean values determined through the Monte 
Carlo simulations, and the errors on each parameter are 
the standard deviation. We then propagate the standard 
deviation on each parameter, along with the correlations 
between parameters, which we find to be much less than 
the standard deviations, to obtain the error bars on the 
far-IR SED model. The degeneracies between the param- 
eters are illustrated in Appendix B. We do find a degree 
of degeneracy between the warm and cold dust tempera- 
tures, but the temperature ranges spanned by each com- 
ponent do not overlap, giving us confidence that we are 
measuring two distinct temperatures. 

We settled on a two temperature model after finding 
very poor fits (high reduced \ 2 values, a factor of ~ 6 
higher) to a single component modified blackbody. The 
introduction of a second dust component produces a sig- 
nificantly better fit, with lower \ 2 values. Since the in- 
dividual photometric uncertainties dominate the noise 
in the far-IR composites, each data point is weighted 
by its associated variance in our fits. In the case of 
the SF SEDs, both sets of data lack photometry from 
~ 18 — 30 /zm. The warm modified blackbody overlaps 
with the mid-IR stacking in the z ~ 1 SF SED, so we 
use this modified blackbody to fill in the gap in the SED 
we were unable to fit. For the z ~ 2 SF SEDs, the warm 
modified blackbody does not overlap with the stacked 
mid-IR spectrum, so we linearly interpolate between the 
stacked photometry and fitted far-IR photometry to fill 
in the area with no photometry. Both of these regions are 
indicated with the dashed lines in Figure[5j We verify the 
far-IR fits by calculating the medians in differential bin 
sizes for the far-IR data using the same techniques de- 
scribed for the near-IR data, and in all cases, our model 
fits match the binned medians. Our composite SEDs for 
each of the four sub-samples are shown in Fig. [5j 

We use longer ground-based (sub)mm measurements 
to verify the accuracy of our composite SEDs. Most of 
our sources are not individually detected in the LABOCA 
and AzTEC+MAMBO maps, and therefore we stack all 
sources in each sub-sample to obtain an average flux 
density at 870 /Ltm and 1.15 mm, respectively. We fol- 



low the (su b)mm stacking method outlined in Pope et 
al. (2008b). The stacking of the SF sub-samples results 
m several > 4cr detections which we plot on the compos- 
ite SEDs in Fig. [5] as the triangles. When the stacking 
results in a non-detection, we plot a 3<r upper limit. The 
(sub)mm detections lie right on the composite SEDs for 
the two SF sub-samples, confirming our modified black- 
body model fits. For the AGN composites, the upper 
limits from stacking are also consistent with our com- 
posite SEDs. 

In Section [4] we discuss the key features of our new 
high redshift SEDs. Given that we report the average 
luminosity, redshift, and stellar mass of the sources that 
go into these templates, they may be useful for fitting 



other high redshift sources, though it is important to 
bear in mind that we created the composites from sources 
that are bright in the mid-IR (5*24 > 200 /iJy). We make 
these new empirical high redshift SED templates publicly 
availabl^ED 

4. FEATURES OF COMPOSITE SEDS 

4.1. Mid-IR spectral features 

We perform mid-IR spectral decomposition on our 
com posi tes as we did for the individual galaxies (see Sec- 
tion |273|). Since the majority of the sources comprising 



both composite SF SEDs are completely dominated by 
star formation, we naturally find that our SF composites 
both have negligible (< 10%) AGN contribution in the 
mid-IR. The featureless AGN composite SED is 100% 
dominated in the mid-IR by an AGN power-law compo- 
nent. The silicate AGN SED is arguably the most inter- 
esting, since it contains both weak PAH features and a 
strong continuum component. Upon decomposition, we 
find the silicate AGN composite to have an 84% AGN 
contribution in the mid-IR. Both the SF SEDs and the 
featureless AGN SED have negligible silicate absorption 
at 9.7 /im, while T9.7 = 0.4 for the silicate AGN compos- 
ite SED. The silicate AGN SED has a power-law slope of 
a = 1.5 from 4 - 15 /im, and the featureless AGN SED 
has a power-law slope of a = 1.6 in the same wavelength 
range. Both numbers are consistent with what is derived 



from local AGN in Mullaney et al. (2011). 



4.2. Dust Temperature and Infrared Luminosity 

Our full composite SEDs are shown together in Fig- 
ure [6] All SED composites peak at approximately the 
same wavelength in the far-IR (~ 100 /im,) which implies 
that each sub-sample has cold dust at approximately 
the same temperature. However, the full SEDs clearly 
show a difference in the relative amounts of warmer dust, 
which we quantify by fitting a two temperature modi- 
fied blackbody model (see Table [2] for the dust temper- 
atures for each sub-sample). A two temperature model 
can be physically understood as follows: the cold dust 
(> 100 /im) is due to the diffuse interstellar medium, 
heated by the underlying stellar population, whereas 
the warm dust originates in SF regions and traces the 
younger stellar population (e.g., Bendo et al. 2012). In 
the case of an AGN, as the AGN becomes luminous the 
dust surrounding it heats up and can radiate in the mid- 
IR and far-IR. Therefore, we might expect to see warmer 
dust in our AGN SEDs than in our SF SEDs. Natu- 
rally, a multi-temperature model is even more ideal, as 
the AGN is not the only heating source for the dust. In 
an AGN still actively forming stars, there should a be a 
third dust temperature that arises from the SF regions. 
Unfortunately, even for our large amount of photometric 
data in the far-IR, fitting more than two dust tempera- 
tures becomes degenerate and produces poor \ 2 values. 
Therefore, we use two temperatures with the assumption 
that these represent the dominant sources of dust radi- 
ation, but there are likely more sources of heating that 
we are unable to separately quantify. 

For the AGN composites, there is a degree of uncer- 
tainty as to where the power-law spectrum ends and the 

21 http:/ /www. astro. umass.edu/~pope/Kirkpatrick2012/ 
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Fig. 5. — Composite SEDs for the mid-IR spectral sub-samples listed in Table ^ In each panel we show photometric (black points) and 
spectroscopic data (black lines) for the sub-sample in the rest frame with our composite SED overplotted in red and the uncertainty of the 
composite shown as the shaded region. The errors on the far-IR data points include both the instrument noise and confusion noise. The 
near-IR and mid-IR data were averaged to obtain the composite while the far-IR data were fit with a two temperature modified blackbody 
model. The blue and green dashed lines illustrate the modified blackbody curves of the individual warm and cold dust components, 
respectively. The open triangles in the upper panels are the stacked (sub)mm detections at observed frame 870 ^m and 1.15 mm. In the 
absence of stacked detections, we plot the 3cr upper limit. In all cases, the (sub)mm detections and upper limits are consistent with our 
composites. 

value, it will produce lower temperatures. We chose to fix 
ft rather than fit it along with the temperatures, but it is 
important to keep the degeneracy of the two parameters 
in mind when interpreting the temperatures. There is ev- 



warm dust component begins. Due to limited spectro- 
scopic data, we are only able to stack spectra below 
~ f 8 /im, and we fit our dust temperatures using only 
far-IR data (above 18 /im). Recent work has suggested 
that the power-law portion of the AGN spectrum c ontin- 



ues to 20 iim, flat tens, and then falls off after 40 /im (Mul 



laney et al. 2011 ). In this case, assigning a single temper- 
ature to this portion of the spectrum is misleading, since 
it is composed of many dust temperatures. However, our 
primary motivation is to compare our AGN composite 
SEDs directly with our SF SEDs, so we choose to fit 
a single blackbody to the warm dust emission around 
^ 20 — 60 /xm. Since the AGN is also radiating into the 
mid-IR power-law portion of the spectrum, the IR lu- 
minosity we derive from the warm dust component only 
will be a lower limit on the AGN emission. 

There is a known degeneracy between the dust emissiv- 
ity, ft, and temperature, such that if we fix ft to a higher 



idence that ft is universal on galaxy- wide scales (Dunne 
fc Eales|200l)j , and as we fix it to the same value for ev- 
ery sub-sample, we do not expect that there is any bias 
created between the sub-samples in this study. 

From our fitting we find that the two AGN composites 
have cold dust temperatures of ~ 33 K, while the SF 
galaxies are only slightly colder at ~ 27 K. However, the 
AGN possess a significantly hotter warm dust component 
(~ 99 K) in comparison to the warm dust component for 
the SF galaxies (~ 57 K). This makes sense if the warm 
dust emission is predominantly heated by the AGN as 
opposed to star formation in dense HII regions. 

To quantify the relative contributions of the hot and 
cold dust components, we calculate Ltr for each compo- 
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Fig. 6. — Comparison of the composite SEDs created from the 
four sub-samples of galaxies. Top panel - No additional normal- 
ization has been applied; the SEDs represent the average intrinsic 
luminosities of sources in each sub-sample. The shaded regions in 
the mid-IR and near-IR correspond to the Icr spread of sources 
around the median luminosities. In the far-IR, the shaded regions 
are computed by propagating the errors on each parameter present 
in the blackbody fitting routine. The z ~ 2 SF SED is the most 
luminous at A > 100 (im, while the featureless AGN SED is the 
weakest at the longest wavelengths. The AGN SEDs are the most 
luminous at 20 fim where the warm dust component dominates the 
spectrum. Middle panel - Composite SEDs have been normalized 
at 100 fim, allowing for more direct comparison of the near and 
mid-IR. No variation of PAH features is seen between SF composite 
SEDs. The featureless AGN SED has the highest mid-IR to far-IR 
ratio. Bottom panel - Composite SEDs have been normalized at 
7 fim. The weak PAH features of the silicate AGN composite are 
visible in direct comparison to the SF SEDs. 



nent separately by integrating over each modified black- 
body model (see Table [3). We compare L^™ and 
^iR d to the total IR luminosity (integrated over 8 to 
1000 /urn), Lir, to determine the fraction of Lir at- 
tributed to each dust component. We then calculate a 
luminosity- weighted average temperature by multiplying 
each luminosity fraction by the respective temperature 
and combining to obtain an effective temperature, T e g 
(Table [1. T eS for the AGN SEDs is more than 20 K 
hotter than T c g for the SF galaxies, demonstrating the 
bolometric significance of the warm dust component in 
AGN-dominated sources. The two SF composites have 
almost the same T e g, indicating little evolution in 7d us t 
with redshift in SF galaxies. 

The Lir values computed by integrating each compos- 
ite SED from 8 to 1000 /um are listed in Table [3) dif- 
ferences between the sub-samples are largely due to our 
flux-limited sample. Lir for the composite SED created 
from z ~ 2 SF galaxies is five times greater than Lir 
for the z ~ 1 SF composite; these luminosities qual- 



ify the z ~ 1 and z ~ 2 SF SEDs to be LIRGs and 
ULIRGs, respectively. The silicate AGN possess a no- 
tably larger Lir (1.7 x 10 12 L ) than the featureless AGN 
(Lir = 8x 10 11 L Q ), though this is likely due to the differ- 
ent redshift distributions of the two sub-samples and not 
any intrinsic evolution of Lir. The luminosities of the 
individual dust components (listed in Table [3]) quantify 
the relative importance of each dust component in the 
different SEDs. For the featureless AGN, the warm dust 
component is more luminous than the cold dust compo- 
nent by a factor of > 2. Though the warm dust does 
not dominate the far-IR SED of any other sub-sample, 
it does contribute at a level of ~ 50% of the total IR 
luminosity. 

This is an important result, since it demonstrates that 
single temperature blackbody models or templates de- 
rived using only longer wavelength SPIRE data will miss 
the im portance of this war m dust component. For ex- 
ample, 



Smith et al 



(2012) derive composite SEDs for 
z < 0.5 galaxies primarily using SPIRE data, which we 
can compare with our composite SEDs. Th e median 
SED of galaxies with L^ ust > 10 11 L from Smith et 
is consistent with our SF composite SEDs in 



al.| ( |2012 

trie mid-IR and beyond the peak in the submm, but it 
is a factor of ~ 2 lower than our SF composites in the 
wavelength range 20 - 60 ^m, where the warm dust dom- 
inates . The discrepancy arises from the fact that [Smith 



et al. (2012) do not possess PACS data for the major- 
ity of their sources and so they may miss a significant 
warm dust component. This comparison underscores the 
importance of full photometric coverage when modeling 
far-IR SEDs. The optimal strategy is to combine deep 
PACS and SPIRE data to obtain a full census of the dust 
emission. 

If, instead of a two temperature modified blackbody, 
we fit our full suite of far-IR data above ~ 20 fim with a 
one temperature modified blackbody, we get a poor \ 2 
fit (reduced x 2 a factor of ~ 6 higher than the reduced 
X 2 values for the 2 temperature fits, due to the fact that 
data > 70 fj.ui is not fit well) , and we calculate an Lir 
that is ~ 30% lower for our SF composites and silicate 
AGN composite. For the featureless AGN composite, us- 
ing a one temperature modified blackbody lowers the Lir 
by ~ 20%. We also derive temperatures that are ~ 20 
K hotter than the cold dust temperatures (see Table [2| 
for all four sub-samples, with the temperatures of the SF 
composites comparable to that of the AGN composites. 
This underscores the importance of the two temperature 
approach, since data points in the range 24 - 100 /im can 
bias any one-temperature fit to a warmer dust tempera- 
ture and lower Lir. 

On the other hand, when photometry spanning the 
range 24 - 100 (im are missing, the dust emission will 
likely be biased to colder temperatures and lower lumi- 
nosities. For our SF composites, only using the cold 
dust for calculations lowers both the Lir and star for- 
mation rates by a factor of ~ 2. The biggest discrepancy 
arises from the featureless AGN, where the warm dust 
is more important than for the other SEDs. There, the 
warm dust accounts for 60% of the total IR luminos- 
ity. When dealing with individual galaxies, often there 
are not enough data points to unambiguously fit both 
the warm and cold dust components. As not account- 
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TABLE 2 

Dust temperatures derived for each sub-sample. 





Effective 


T 


T 


Sub-sample 


Temperature a 


(warm dust) b 


(cold dust) b 


2 ~ .1 SF galaxies 


37 ± 3 K 


55 ± 6 K 


25 ± 2 K 


z ~ 2 SF galaxies 


40 ± 3 K 


59 ± 5 K 


28 ± 2 K 


Silicate AGN 


58 ± 2 K 


98 ± 2 K 


35 ± 3 K 


Featureless AGN 


65 ± 2 K 


100 ± 1 K 


33 ± 1 K 



Note. — See Table [3] for IR luminosities. 
a Calculated from a luminosity-weighted average of the warm and cold 
dust components. See text for details. 

b Temperatures are the mean values determined from Monte Carlo 
simulations (see fj3), with the errors being the respective standard 
deviations 



ing for any warm dust may significantly bias results, it 
is perhaps best to use templates which account for both 
important dust components as observed in high redshift 
galaxies. 



4.3. Overall SED shape 

In the top panel of Fig. [6j which compares the intrinsic 
luminosities of the sources going into each composite, 
the AGN SEDs have more emission than the SF SEDs 
in the range 15 to 30 /tin. At z ~ 2, the MIPS 70 /im 
passband covers this wavelength range. In fact, for our 
high redshift sources (z > 1.5), 64% of those detected 
at 70 /im are AGN dominated. For our lower redshift 
sources (z < 1.5), only 19% of 70 /im-detected sources 
are dominated by an AGN in the mid-IR. At the peak 
of the SED, the most luminous sources are the z ~ 2 SF 
galaxies and silicate AGN, reflecting their higher total 
IR luminosities listed Table |U 

We normalize the composite SEDs to 100 and 7 /im in 
the middle and bottom panels of Fig. [6j respectively, 
to allow a more direct comparison of the mid-IR spec- 
tral features. The mid-IR slopes of the AGN SEDs are 
remarkably similar, despite the absence/presence of fea- 
tures, but the featureless AGN have a much larger mid- 
IR to far-IR ratio. In fact, the featureless AGN SED ef- 
fectively flattens out at A > 20 /im. Both AGN SEDs lack 
a visible stellar bump at 1.6 /im, indicating that the rest- 
frame near-IR emission is also dominated by the AGN. 
The SF and silicate AGN SEDs all exhibit PAH features, 
though they are much weaker for the silicate AGN SED. 

In the middle and bottom panels of Fig. |6j we can 
directly compare our two SF composite SEDs to look for 
any evolution in SED shape with redshift and/or Lir. 
The PAH features of the two SF SEDs are nearly iden- 
tical in strength and shape. There appears to be no 
evolution of PAH features between LIRGs at z ~ 1 and 
ULIRGs at z ~ 2. The far-IR part of the SF SEDs are 
also remarkably similar in shape, as quantified in the av- 
erage dust temperatures (Table pi. The stellar bump at 
1.6 /im is more prominent relative to the IR emission for 
the z ~ 1 SF SED than for the z ~ 2 SF SED, reflecting 
the higher average specific star formation rate of galaxies 
at z ~ 2 compared to z ~ 1, which we discuss in the next 
section. 



4.4. Star Formation Rates 



We calculate the star formation rate (SFR) for each 
composite SED according to the formula 



SFR 

M Q yr- 



1.72 x 10 



-10 



rSF 



(2) 



which assumes a Salpet er IMF (Salpeter | 1955[ ) and con- 
tinuous star formation ( Kennicutt||l998| )" INote that we 
only want to use the portion of .Lir that is heated by 
SF activity when using this equation, otherwise we will 
over-predict the SFR. For the SF composite SEDs, we 
assume that both the cold and the warm (~ 60 K) dust 
components are heated by star formation, since there is 
no evidence of AGN activity in the composites based on 
the mid-IR spectrum. 

We calculate the fraction of the AGN SED that can be 
accounted for by the SF SED template when both are 
normalized at 200 /im, beyond which all the dust emis- 
sion is presumed to be due to SF. We calculate the Lir 
of the normalized SEDs and take the ratio of Lir for the 
z ~ 1 SF SED and featureless AGN SED, and for the 
z ~ 2 SF SED and silicate AGN SED. We use the SF 
composite SED that most closely matches the total Lir 
of the respective AGN composite when determining the 
Lir ratios. We find that SF activity accounts for 56% of 
the Lir for the silicate AGN SED and 21% of the feature- 
less AGN SED. We therefore scale the L m by 56% and 
21% for the silicate and featureless AGN, respectively, 
when calculating the SFR. 

Unsurprisingly, our z ~ 2 SF composite SED has the 
highest SFR: 344 M Q yr" 1 . The z ~ 1 SF composite SED 
has an SFR of 73 M yr" 1 . The silicate AGN SED has 
an average SFR of 159 M yr -1 , whereas the featureless 
AGN SED has a much lower average SFR of 28 M Q yr" 1 . 
If we instead assume that all of the Lir in the AGN 
composites is heated by SF, we would over-predict the 
SFR by a factor of - 2 - 4. 

As a consistency check, for the two AGN SEDs, we con- 
servatively attribute the warm dust at 100 K to the AGN 
and the cold dust to the SF activity and recalculate the 
SFR using Ljg ld instead of the scaled Lir; naturally, this 
conservative estimate should provide a lower limit on the 
actual SFR, since there is likely some portion of warmer 
dust emission that arises from SF. Using just the cold 
dust component, we find the same SFR of 28M Q yr _1 
for the featureless AGN, and an SFR of 141 M Q yr -1 for 
the silicate AGN, which are very similar to the SFRs we 
calculated above, so most of the SF activity in each AGN 
composite is accounted for by the cold dust component. 
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There is a tight correlation betw een stellar mass and 



SFR which evolves with redshift ( Daddi et al.||2007 
Dunne et all [20091 [Pannella et a"L"1|2009a[ |Elbaz et al., 
2011[|Lin etaL| 2012), and as such, an interesting compar- 
ison is given by the specific star formation rates (sSFR), 
defined as: sSFR = S FR/M^. We list the median stellar 
mass (see Section |2.5[ ) and sSFR of each composite SED 
in Table [3J For the SF composites, which have compara- 
ble median stellar masses, we see an increase in the sSFR 
from z ~ 1 to z ~ 2, consistent with w hat is observed for 
all SF galaxies (e.g., |Elbaz et al.|2011[ ). The silicate AGN 
and z ~ 2 SF SEDs have comparable stellar masses and 
redshifts, but the sSFR for the z ~ 2 SF SED is more 
than a factor of 2 higher. Similarly, the featureless AGN 
composite has a higher stellar mass and ~ 4 times lower 
sSFR than the z ~ 1 SF composite. As we will discuss in 
Section 5] these values are consistent with a scenario in 
which the AGN sources are in a later phase in the evo- 
lution of IR luminous galaxies, after the active period of 
star formation has been quenched. 

4.5. Comparison of New High z SEDs to Local 
Templates 

Prior to the era of Herschel, far-IR (40 - 200 /mi) data 
sampling the peak of the SED were rare, particularly for 
high redshift sources. Estimates of bolomctric luminosity 
were performed by fitting locally derived templates to 
mid-IR and/or submm data and extrapolating to the far- 
IR. With our wealth of far-IR data from Herschel, we are 
in a position to compare the commonly applied locally 
derived templates with our new empirical high-redshift 
SEDs. 

In Figure [7] we compare our SF composite SEDs to 
the SEDs of local galaxies often used as standard tem- 
plates in the literature. These local templates come from 
combining all known data on known sources including 
available IRAS (12, 25, 60, 100 /im), Spitzer (24, 70, 
160 /im), and SCUBA (85 Atm) photometry as well as 
mid-IR spectrosco py (e.g., Forster Schreiber et al.|[2003 
Armus et al.|2007 1 . The local templates lack photometry 
trom 160 - 850 /tin, making it difficult to constrain the 
cold dust emission. 

In the top panel of Fig. [7j we compare our z *~ 1 
composite SF SED to the SED of local starburst M 82. 
Though M 82 has a lower L m (~4x 10 10 L Q ) than our 
z ru 1 SED, it is worthwhile to compare the two SEDs, 
since our mid-IR classification is based on the M 82 tem- 
plate. We normalize M 82 to our composite SED at 
7.0 /im. Our z ~ 1 composite SED almost exactly re- 
produces the shape and strength of the PAH features of 
M 82, affirming the reliability of comparing the mid-IR 
PAH features of this local galaxy with the mid-IR spectra 
of high redshift sources. However, M 82 has less far-IR 
emission, and it is weighted towards warmer dust than 
our z ~ 1 SF composite, indicating that any mid-IR sim- 
ilarities do not carry over into the far- IR. 

The bottom panel of Fig. [7] shows our z ~ 2 composite 
SF SED compared to two local prototypical (U)LIRGs, 
Arp 220 (L m = 1.4 x 1O 12 L ) and NGC 6240 (L m = 
6.3 x 1O 11 L ), both normalized to our composite SED 
at 7.0 /im. While the mid-IR spectrum of NGC 6240 is 
very close to that of our high redshift composite, Arp 220 
shows substantially more silicate absorption. The promi- 
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Fig. 7. — Comparison of our sources and composites with local 
SEDs, overplotted on the photometric (gray points) and spectro- 
scopic (gray lines) data for a given SF sample. Top - z ~ 1 SF com- 
posite SED compared to the SED of M 82 (normalized at 7.0/im). 
M 82 is in excellent agreement with our z ~ 1 composite SED 
in the mid-IR, but is weighted towards warmer dust temperatures 
and underpredicts the luminosity in the far-IR/submm. Bottom 
- z ~ 2 SF composite SED compared to the SEDs of Arp 220 
and NGC 6240 (both normalized at 7.0 fim). The local and high z 
ULIRG SEDs show fairly good agreement in the range 5-11 jim, 
but both local SEDs are too luminous in the far-IR, particularly at 
the shortest wavelengths, indicating more warm dust. In addition, 
Arp 220 has stronger silicate absorption than our composite SED. 



nent absorption feature at 9.7 /im is possibly attributabl e 
to AGN activity in Arp 220 (e.g. 
while our z 



Armus et aL||2007 ) 
2 composite SF SED has a negligible AGN 
contribution. Neither of these local (U)LIRGs match our 
z ~ 2 SED in the near or far-IR; specifically, Arp 220 and 
NGC 6240 peak at shorter wavelengths in the far-IR, in- 
dicating much more warm dust emission. It is possible 
the difference in warm dust emission arises from the fact 
that local ULIRGs are merger-induced, triggering AGN 
activity, while our sampl e of high redshift ULIRGS may 
not be (see Section 5.2.1 ). These local ULIRG templates 
differ significantly from our z ~ 2 ULIRG composite, and 
in fact, there are not even any individual sources that go 
into the composite that fit the local ULIRGs from 20 
- 80 /im (see the gray points showing all galaxies that 
go into the composite). This illustrates the difficulty of 
applying such local templates to high redshift systems. 

We compare our silicate AGN composite SED with the 
SEDs of local AGN Mrk 231 (L m = 3.2 x 10 12 L Q ) and 
NGC 1068 (Lm = 1-6 x 10 U L Q ) in Figure § Both 
Mrk 231 and NGC 1068 have been normalizca to our 
composite at 7.0 /im. Our composite and Mrk 231 ap- 
pear to have a similar 9.7 /tm silicate absorption feature 
and a similar slope in the mid-IR, but Mrk 231 peaks 
at shorter far-IR wavelengths than our composite SED, 
indicating more warm dust. NGC 1068 looks fairly sim- 
ilar in shape in the far-IR to our composite, though it 
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TABLE 3 

Luminosities and other parameters derived for each sub-sample. 



Sub-sample 


Ltr (8-1000 (im) 
(10 12 L Q ) 


r warm a 

(10 11 L ) 


r cold a 
*TR 

(10 11 Lq) 


r cold 
%R 
/ total 
■°IR 


(10 11 L ) 


T total 
^IR 

L 8 


M, 
(10 11 M ) 


SFR 

(M yr- 1 ) 


sSFR 
(Gyr- 1 ) 


z ~ 1 SF galaxies 


0.42 ± 0.17 


2.0 ± 0.5 


2.0 ± 0.6 


0.46 


0.65 


6.57 


0.74 


73 ± 29 


0.99 


z ~ 2 SF galaxies 


2.00 ± 0.71 


8.7 ± 1.9 


10.2 ± 2.1 


0.51 


2.53 


7.92 


1.17 


344 ± 122 


2.94 


Silicate AGN 


1.65 ± 0.54 


7.0 ± 0.5 


8.5 ± 2.4 


0.50 


2.60 


6.35 


1.15 


159 ± 52 c 


1.38 


Featureless AGN 


0.76 ± 0.07 


4.4 ± 0.2 


1.6 ± 0.3 


0.21 


1.94 


3.94 


1.23 


28 ± 3 C 


0.23 



a Determined by integrating under the modified blackbody curve for the warm and cold components separately. 

b Measured by integrating under the IRAC 8 /im filter for each composite SED. 

c Calculated using a scaled for the AGN SEDs (shown in italics to distinguish). 
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Fig. 8. — Comparison of our silicate AGN composite SED with 
SED templates for the local AGN Mrk 231 and NGC 1068. The 
templates are overplotted on the photometric (gray points) and 
spectroscopic (gray lines) data for the silicate AGN sample. Both 
Mrk 231 and NGC 1068 have been normalized to the silicate AGN 
SED at 7.0 /im. Mrk 231 peaks at shorter wavelengths in the far-IR 
than our template, while NGC 1068 falls off sharply below 3 fim. 

is intrinsically less luminous. Furthermore, NGC 1068, 
a local Compton-thick AGN, has less silicate absorption 
in the mid-IR and falls off sharply below 3 /im, differing 
from both our high redshift SED and Mrk 231. 

With the emerging abundance of Herschel photometry, 
recent studies have begun to demonstrate the unsuitabil- 
ity of locally-derived templates in fitting high redshift 
alaxies using longer wavelength PACS and SPIRE data 



Dannerbauer et al 



on et al. ||2ljl0|). 



20101 |Elbaz et al.l|2010| |20Tlj \Nor- 
hese p revious results fit templates 
from Cfiary & Elbaz ( 2001| CE01) templates to individ- 



ual galaxies with far-IR photometric data points, which 
each have significant uncertainty. In this study, we are 
able to compare the locally-derived CE01 templates to 
templates derived from high-redshift data for a statis- 
tical sample of galaxies. We plot the CE01 templates 
with the Ltr corresponding to each of our SF compos- 
ite SEDs in Figure [9j with no renormalization. For the 
z ~ 1 SF composite (top panel), the CE01 template 
reproduces fairly well the near and mid-IR portions of 
the composite SED. However, it peaks at slightly shorter 
far-IR wavelengths and has less cold dust emitting at 
(sub)mm wavelengths. For the z ~ 2 SF composite 
(bottom panel), the CE01 template fails to reproduce 
any portion of our z ~ 2 composite SED. The mid-IR 
spectrum of the CE01 template has a strong continuum 
component and weaker PAH features than we observe in 
the average z ~ 2 ULIRG of similar luminosity. Again 
the far-IR emission peaks at shorter wavelengths, corre- 
sponding to warmer dust emission, than our z ~ 2 SF 
SED. The stacked (sub)mm points (plotted as the tri- 
angles) are also inconsistent with the CE01 templates, 
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Fig. 9. — Comparison of our z ~ 1 (top) and z ~ 2 (bottom) 
SF composite SEDs with the locally-derived CE01 templates of the 
same luminosity. The templates are overplotted on the photometric 
(gray points) and spectroscopic (gray lines) data for a given SF 
sample. For the z ~ 1 SF galaxies, the CE01 template traces the 
composite SED in the near and mid-IR but differs in the far-IR. 
For the z ~ 2 SF galaxies, the CE01 template fails to reproduce 
any portion of the composite SED. 

confirming that our high redshift sample has more cold 
dust emission than the CE01 templates. 

Overall, we find that most local templates fail to accu- 
rately reproduce both the dust temperatures and mid-IR 
spectral features of our high redshift composite SEDs. 
It is possible that things will be more consistent once 
Herschel data are incorporated to create more complete 
SEDs of local galaxies and better constrain the cold dust 
emission. Otherwise, this is evidence that the SEDs of 
IR luminous galaxies evolve with redshift. The disparity 
of the locally derived templates from our high redshift 
composites illustrates the need for caution when apply- 
ing local templates to high redshift galaxies. 

5. DISCUSSION 

In Sections [3] and |4j we presented our high redshift 
composite SEDs and outlined key differences and simi- 
larities between the SF- and AGN-dominated IR lumi- 
nous galaxies. We also compared our new high redshift 
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SED templates to those of local galaxies of similar lumi- 
nosities and found them to be significantly different. We 
now discuss what these results might tell us about the 
evolution of IR luminous galaxies. 

5.1. Quantifying the AGN Component 

We begin by exploring the X-ray and mid-IR properties 
of the AGN comprising our composite SEDs. Though an 
in-depth X-ray analysis of our AGN sample is beyond 
the scope of this paper, we briefly discuss two important 
X-ray parameters commonly used to interpret AGN: X- 
ray detection fraction and X-ray luminosity. Our entire 
sample of AGN sources (38) have a 58% X-ray detection 
fraction. When we classify our sources into silicate and 
featureless AGN, only 9/22 (41%) silicate AGN are in- 
dividually detected in the deep Chandra X-ray imaging, 
whereas 11/12 (92%) featureless AGN are detected in 
the X-ray. The four unclassifiable AGN have a detection 
fraction of 50%. 

We calculate the in trinsic X-ray luminosit y, L2-iokcV, 
using equation (4) of Mullaney et al. (2011). 



log 



{ 10 43 



t AGN 
^IR 



ergs 



= 0.53 + 1.11 log 



/ ^2-lOkeV 

V10 43 ergs" 1 



(3) 

To disentangle the luminosity due to the AGN from the 
host galaxy, we follo w th e same method used to calcu- 
lated the SFR (see \ AA). We normalize all composite 
SEDs to 200 /im and calculate the ratio of the normal- 
ized L m for the z ~ 1 SF SED and featureless AGN SED, 
and for the z ~ 2 SF SED and silicate AGN SED. We find 
that AGN activity accounts for 44% of the Ltr for the 
silicate AGN SED and 79% of the featureless AGN SED, 
and accordingly we scale each respective Lm by these 
amounts when calculating £ 2 _iokcV- This gives an in- 
trinsic X-ray luminosity of £2-10 kcV = 5.30 x 10 44 erg s _1 
for the silicate AGN and L 2 -iokcv = 4.46 x 10 44 ergs _1 
for the featureless AGN. 

We look for evolution between our high redshift AGN 
and local AGN by comparing to the sample of moder- 



at e luminosity (La-mimy 



10 4 ^ - lO^ergs- 1 ) AGN 
B y cros s-matching the Swift- 



Mullaney et al.| J201 
BAT sample ([Tueller et"aLlT2008t of X-ray~A GN with 



the Spitzer-IRS archive, Mullaney et al. (2011) identified 
a sample of nearby (i.e., z < 0.1) galaxies whose mid- 
infrared spectra are strongly AGN-dominated. By care- 
fully subtracting any contribution from the host galaxy 
at far-IR wavelengths, they were able to empirically de- 
fine the infrared emission due solely to the AGN in these 
local galaxies and derive an average, intrinsic AGN SED 
at 6 - 100 ^m. In Figure[l0j we plot the mean AGN SED 
and the high and low luminosity AGN SEDs from the lo- 
cal population along with our high redshift AGN compos- 
ite SEDs. The high and low luminosity AGN SEDs were 
derived by averaging the SEDs of log(L 2 _iokcv) > 42.9 
and log(Z/2-iokcv) < 42.9 AGN, respectively. The low 
luminosity AGN template and the mean AGN template 
turn over around 40 — 50 fiin whereas the high luminosity 
AGN template turns over at 32.0 /Ltm. In both of our high 
z AGN composites, the warm dust component, which we 
attribute to the AGN, turns over at approximately the 
same wavelength (29 /im) as the high luminosity local 
AGN template. 
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Fig. 10. — We compare our compos ite AGN SEDs (solid lines) 
with the average AGN SEDs derived in Mullaney et al. (201ll, nor- 
malized at 7 fim. The blue dashed line shows the mean AGIN SED, 
while the green and orange dot-dashed lines show the high and 
low luminosity AGN SEDs, respectively, which have been created 
using the SEDs of local AGN. We overplot the warm dust compo- 
nent for each of our composite SEDs as the dotted line. The warm 
dust emission peaks at the same wavelength as the high luminosity 
AGN SEDs, which is consistent with the high L^— iokev values we 
derive for our composites. 



The mid-IR spectral decomposition of our AGN com- 
posite SEDs resulted in an 84% contribution of the AGN 
to the mid-IR luminosity for the silicate AGN SED and 
a 10 0% contribution for the featureless AGN SED (see 
£4.1). If we make the simplifying assumption that the 



warm temperature component is due entirely to the AGN 
and the cold component to the host galaxy, then we can 
quantify the contribution of the AGN to the total IR lu- 
minosity by subtracting the contribution of Lf^ from 
L m (see Table [|). For the silicate AGN SED, we find 
an IR contribution of 50%, and for the featureless AGN, 
a contribution of 79%, in good agreement with what we 
find when we remove the host component by scaling the 
SF composite SEDs. Both SEDs exhibit a lower overall 
contribution from the AGN component to Lib. than that 
indicated by the mid-IR alone, since the SF component 
dominates at the longest IR wavelengths. As a check on 
the validity of our simple assumptions about the two tem- 
perature compo nents, we also compare the in^* N that we 
derive from the Mullaney et al.| (12011 ) mean AGN tem- 
plates (Fig. [TO]). Fitting the lMullaney et al.|(|2011[ ) mean 
AGN templates to the mid-IR of our AGN SEDs, we find 
an AGN contribution to Ltr, of 62% and 83% for the sil- 
icate and featureless AGN, respectively. These results 
are remarkably consistent with our modified blackbody 
fits, which is reassuring. With our blackbody fits, we are 
able to attach physical parameters, namely dust tem- 
peratures, to each population which eases comparisons 
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between our samples and with other galaxy populations 
at low and high redshift. 

5.2. Evolution of infrared luminous galaxies 

In the local Universe, all ULIRGs and most LIRGs 
are observed to be caught in the act of a major merger 
(Sanders & Mirabel 1996), where the intera ction is the 



al. 


1996 


Bushouse et al. 


2002) 



A popular scenario for 
linking IK luminous g alaxies dominated by SF and AGN 



activity is outlined in Sanders et al. (119881. During the 



early stages of a merger of two disk galaxies a dusty star 
burst is triggered, but as the black holes merge together 
and are fed by the disrupted gas, an AGN begins to domi- 
nate the IR emission and destroys or blows away the dust 
and gas, quenching the star formation. Eventually the 
AGN runs out of fuel, and the galaxy settles down as 
a massive elliptical galaxy. This scenario is a plausible 
explanation for local ULIRGs. 

At z ~ 1-3, the situation is different. While some 
high red shift ULIRGs are major mergers (e.g 



has a lower sSFR than the z ~ 2 SF composite SED. 
The AGN properties of our silicate AGN SED are con- 
sistent with a later evolutionary stage than our z ~ 2 SF 
galaxies, regardless of whether the evolution is driven by 
major mergers or some other process. 

5.2.1. Main sequence of star forming galaxies 

One way that has been proposed to determine if a high 
redshift galaxy is undergoing a starburst event triggered 
by a merger is using the ratio of SFR to stellar mass. 
Several recent studies have found a tight correlation be- 
tween SFR and stellar mass which has been called the 
'main sequence' for star-forming galaxies (e.g 



Noeske 



et~aL|[2007l JElbaz et al.|[2007[ |Daddi et al.l|2007| [Magdis 
et al.||2 0To| . Furthermore, this main sequence shifts to 
SFR from z 



to an increase in the 
If 



2011) 



al. 



Engcl et 



2010), many high redshift ULIRGs show no evidence 



for any merger or interactio n, from studies of their mor- 



phology or dynamics (e.g., Daddi et al. 2010 Tacconi 



et al. 



20101. The higher gas tractions and longer gas (20111 and Eqn. 2 in 



higher SFR from z = - 3, leading 

average sSFR with redshift (e.g., [Elbaz et al i 
galaxies on the main sequence are undergoing continuous 
star formation then starburst galaxies will lie above the 
main sequence, their SFRs possibly boosted by a merger 
or interaction. 

Both of our SF composites and our silicate AGN com- 
posite have sSFRs that place them on the sS FR-redshift 
main s equence defined by Eqns. 13 an d 24 in l Elbaz et al.| 



consumption timescales suggest that the intense star for- 
mation in these galaxies is fuel ed continuously, p erhaps 
through cold gas strea ms (e.g. 



Dekel et al. 2009b. Sev- 



eral recent studies (e.g., |Schawinski et al.|20il[|Mullaney| 
et al.||2012| have found that many moderate luminosity 
fo 



Pannella et al | ( |2009ap . The fea- 
tureless AGN, however, lie decidedly below both of these 
relations, though not in the regime designated as "qui- 
escent" (< 10% of th e average sSFR at a given redshift, 
Mullaney et al.|[2012 ) . Our stellar masses are a factor of 



2 lower than those in 



AGIN are found in high redshift IR luminous galaxies 
that are not necessarily undergoing a major merger, sug- 
gesting that internal instabilities could be the primary 
mechanism fueling the AGN. 

We found significant differences in the IR SED between 
our mid-IR identified SF and AGN galaxies; specifically, 
the AGN sources have effective dust temperatures ~ 20 K 
higher, indicating that the luminous AGN is not only re- 
sponsible for heating the hot dust in the mid-IR, but also 
has an impact on the warmer far-IR dust and can even 
dominate the bolometric luminosity (e.g. our featu reless 
AGN in Table [3j see also |Di'az-Santos et al."1|2010[ ). It is 
tempting to try to place our SI 1 '- and AGN-doniinated 
galaxies into an evolutionary sequence, where the SF 
galaxies represent an early phase of active star forma- 
tion, and the AGN galaxies represent a later phase once 
the black hole has been fueled and is able to heat the 
dust to warmer temperatures and quench the star for- 
mation. This evolutionary sequence may or may not be 
triggered by a major merger at high redshift. Any direct 
evolutionary comparison can only be tentatively drawn 
between our z ~ 2 SF SED and our silicate AGN SED, 
since these two SEDs are comparable in Lir and stellar 
mass, and both are primarily composed of galaxies with 
a redshift of ~ 2. There is some tentative evidence that 
a major merger prod uces more warm dust emission than 
normal SF activity (Hayward et al. 2012). Our z ~ 2 
SF SED exhibits much less warm dust emission than the 
local ULIRG SEDs, indicating high redshift ULIRGs are 
not analogs to local ULIRGs and may follow a different 
evolutionary path in which they are normal SF galaxies 
until the AGN at the center begins to dominate the bolo- 
metric luminosity. Our silicate AGN composite SED is 
luminous in the X-ray (£2-iokcV ~ 5 x 10 44 ergs _1 ) and 



Elbaz et a l. (2011) due to differ- 



ent stellar population models (see J 275 ) , so if we increase 



our stellar masses by this amount, our sSFRs will conse- 
quently be a factor of 2 lower, placing them even further 
a way from the SB re gime . 



Elbaz et al. 



( 2011[) discovered another tight correlation 



between the 8 /jm luminosity, Lgum; and Lir, mimick- 
ing the SFR-Af* main sequence, and defined an alter- 
nate main sequence using the parameter IR8 = Xir/L 8 . 
IR8 is calibrated to distinguish between starbursting and 
main sequence SF galaxies, but it has been shown to be 
useful when applied to AGN as well. We calculate the 
8 /Urn luminosity of our SF and AGN sources by applying 
the IRAC 8 fim passband to each composite SED in the 
rest frame. We calculate IR8 = Ljr/Ls for each com- 
posite and list the values in Table [3| We do not distin- 
guish between main sequence and starburst galaxies for 
individual sources since we are interested in the average 
properties of our sub-sam ples. Using a large sample of 
GOODS- Herschel galaxies, |Elbaz et al.l ( |2011| ) define the 
main sequence to be IR8 = 2.7 to 7.8 (range is the 68% 
dispersion around the mean IR8 value of 4.9). Interest- 
ingly, IR8 for both of our AGN composites fits comfort- 
ably within this range. The z ~ 1 SF SED has an IR8 
of 6.6, which is well within the main sequence, while the 
z ~ 2 SF SED, with an IR8 of 7.9, lies just above the 
main sequence range. Although most of the galaxies in 
our z ~ 1 SF sub-sample are main sequence galaxies, the 
z <~ 2 SF sub-sample may contain a mix of both main 
sequence galaxies and starbursts, according to the IR8 
parameter. 

During the active stages of a major merger, the AGN 
produced is likely to be heavily obscured. As our present 
study does not rely on X-ray detection to identify AGN, 
we are not biased towards unobscured AGN. Therefore, if 
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Fig. 11. — Our z ~ 1 and z ~ 2 SF SEDs (top and bottom panel, 
respectively) compa red against the ma in sequence and starburst 
templates derived in |Elbaz et al.| ( |2011[ l . Both the main sequence 
and starburst templates overpredict en 
wavelengths (> 200 /im). 



emission at far-IR/submm 



AGN at high redshift primarily reside in (U)LIRGs trig- 
gered by a major merger, our main-sequence/starburst 
classification should reflect this. As neither the silicate 
AGN nor featureless AGN lie in the starburst regime, 
according to the IR8 criterion and the sSFR-z relation, 
our findings are in agreement with recent studies that 
suggest the majority of AGN at high redshift (z ~ 1 - 3) 
reside in normal, main sequence galaxies (e.g. , Lutz et al. 

However, 



2012 Xue et al. 2010). 



2010| |Mullaney et al 

two caveats must be taken into consideration. First, a 
compact starburst and an AGN can, in effect, wash each 
other out in the IR8 parameter. The starburst can in- 
crease the ratio of the far-IR to the mid-IR, but the AGN 
will increase the amount of mid-IR luminosity. Together, 
these effects act to produce a main-sequence IR8 value. 
Second, both main sequence criteria may be limited in 
usefulness when discussing AGN. A major merger may 
immediately trigger a starburst, but may take more time 
to funnel material to the AGN, and this can produce a 
delay between starburst activity and a luminous AGN on 
the order of ~ 10 s yrs. It is possible that any signatures 
of a merger have disappeared before the galaxy displays 
strong AGN signatures. 

Based on the IR8 parameter, most of our AGN and 
SF galaxies do not show enhanced IR luminosity and 
may not be a phase of the popular local ULIRG major 
merger scenario, although they might still be linked in 
some other evolutionary sequence (e.g. involving minor 
mergers or cold gas accretion). Naturally, our findings 
are only for the average of all of our high redshift sub- 
samples of galaxies and do not necessarily apply on an 
individual basis. 



Elbaz et al. (2011) present templates for typical main 



ing photometry for main sequence and starburst galaxies 
at z < 2.5 from 3 - 350 /im. We compare both of our SF 
SEDs to the main sequence and starburst templates in 
Figure |11| Upon comparison to ou r SF templates, we 
find tha/tooth the Elbaz et al. (2011 ) main sequence and 
starburst templates overpredict the amount of emission 
at wavelengths greater than 200 /jm rest-frame exhibited 
by both our z ~ 1 and z ~ 2 SF templates, indicat- 
ing even more cold dust emission. This is possibly due 
to a difference in the average intrinsic IR luminosity or 
redshifts for each of these templates. In addition, there 
are differences in the way the templates were created, 
which could also account for some of the variation in cold 



dust emission. [Elbaz et al. (1201 lh only include the long 



wavelength SPIRE photometry m their fit when there 
is a detection in the prior catalog whereas we include a 
SPIRE measurement for every galaxy regar dless of if it 
is a formal detection (see fj3|. Furthermore, Elbaz et al.| 
( 2011 1 fit the far-IR data for individual sources (requiring 
at least one measurement above 30 rest wavelength) 
with the CE01 templates to calculate each Ltr, and then 
normalized every source to the same Ltr before stacking 
to create their templates. This can introduce a degree 
of uncertainty since it requires an extrapolation between 
the data points. We normalize out sources in the mid-IR 
where we have the IRS spectra so we don't need to rely 
on any templates. The amplitude and shape of our SEDs 
at these longer wavelengths is validated by our stacked 

_ (|20njrtemplates 
ikely due to the fact that we 



submm measurements (see triangles in Fig. |11|). Most of 
the mismatch between the |Elbaz et al 
and ours in the submm is 



are com paring sources at diff erent redshifts and luminosi- 
ties; the Elbaz et al. (20111 SEDs include all galaxies at 
z < 2.5 and will be weighted towards galaxies at z < 1 
with lower average Lm,, which may indeed be intrinsically 
cooler, than our SF composite SEDs. As shown earlier in 

> 1 galaxies, 



Section 4.5 



local templates do not fit the z 
and we expect the SEDs of star forming galaxies to start 
to evolve somewhere between z ~ 1 and z ~ 0. A more 
in depth study of the SEDs of star forming galaxies at 
z < 1 is needed to determine when the SED shapes begin 
to differ from local templates. 

5.2.2. Obscuration in AGN 

Our AGN sources are clearly split into two sub- 
samples, those with a measurable silicate feature and 
those without, and the same is seen in samples of local 
AGN. Silicate absorption has been linked wi th the most 
heavily obscured, Compton-thick AGN (e.g 
2011 Georgantopoulos et al. 2011 ) 



Sales et al. 
Sources with absorp- 
tion have an optically thick, smooth dust distribution, 
wher eas clumpy dust produces only shallow absorption 



(e.g.,|Spoon et al.||2007l iNenkova et al.||2008), although 



recently Goulding et al.| ( |2TFrZ ) lind that the deep silicate 



sequence and star bursting galaxies, computed by stack- 



absorption features seen in local Compton-thick AGNs 
are predominantly ca used by extinction due to the host 
galaxy. According to Levenson et al. (2007), based on 
modeling local AGN, the strength of the silicate absorp- 
tion feature in our silicate AGN SED does not neces- 
sarily require the SED to be heavily obscured; neverthe- 
less, there are significant differences between the silicate 
AGN and featureless AGN subsamples. Interestingly, in 
our sample only 41% of silicate AGN are individually 
detected in the deep Chandra X-ray imaging, whereas 
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92% of featureless AGN are detected in the X-ray. This 
supports the idea that the silicate AGN are more heavily 



obscu red, possibly even Compton-thick (lAlexander et al 
20081, on an individual basis, while the featureless AGIN 
are less obscured and can have a stronger effect on the 
host galaxy. A more detailed study of the individual X- 
ray properties of our IRS sample is reserved for a future 
paper. 

Broadly speaking, silicate AGN may represent the 
transition stage from a star formation do minated galaxy 
to an unobscured AGN, following the |Sanders et al.| 
(1988) prescription for creating luminous QSOs. For our 
AGJN sources, we assume the mid-IR emission arises from 
the hot dust torus surrounding the AGN, while the emis- 
sion at > 60 /im arises from the cold dust heated by star 
formation in the host galaxy. The persistence of signifi- 
cant cold dust observed in our silicate AGN suggests that 
these nuclei are not yet having a strong effect on the star 
formation occurring in the host galaxy. For an unob- 
scured AGN, the mid IR emission from the warm dust in 
the torus will dilute the dust emission in the surround- 
ing galaxy, so the estimated amount of cold dust is sup- 
pressed or overrun by the warm dust. Currently, a precise 
picture of AGN feedback has not been proven, but one 
prominent line of thought is that as the AGN becomes 
more powerful, the feedback from the AGN is responsi- 
ble for actively shutting down star formation. AGN feed- 
back could be responsible for the lack of any SF features 
or silicate absorption in the mid-IR, the suppression of 
cold dust emission, and the much lower SFR observed 
in our featureless AGN compared to our silicate AGN. 
Alternatively, the featureless AGN sources might just be 
more luminous AGN, such that they outshine any of the 
host galaxy emission, although our infe rred X-ray lumi- 
nosities are inconsistent with this (see [5.1 1. Of course, 
while these scenarios are possible in general, based on 
the average properties of our sub-samples of galaxies, it 
is impossible to state that this is the case for every galaxy 
within each sub-sample, since we are probing a range of 
redshifts and Lir. Furthermore, with the level of silicate 
absorption we see in our templates, it is difficult to dis- 
tinguish between a clumpy and smooth torus geometry. 
In this case, it is possible that we are simply seeing iden- 
tical AGN at different viewing angles which produces the 
different shapes in the spectra. 

While our SEDs are consistent with a scenario in which 
the featureless AGN are in a later evolutionary stage than 
the silicate AGN, and are having more of an impact on 
the host galaxy, evidenced by the lack of mid-IR fea- 
tures, the suppression of cold dust, and the higher X-ray 
detection fraction, we cannot rule out other explanations 
for the differences between the SEDs. Detailed modeling 
of emission from different AGN components that con- 
tribute to our composite SEDs is beyond the scope of 
this paper. Fully disentangling why IR luminous AGN 
at high rcdshift have noticeably different spectral shapes 
is a necessary step to producing a complete picture of 
galaxy evolution. 

5.3. A closer look at our high redshift SEDs 

Looking only at the high redshift SF sources, we might 
ask whether there is evidence for evolution in the SEDs 
between z ~ 1 and z ~ 2. We showed in Fig. [6] that 
when normalized in the mid-IR or far-IR, our two SF 



composite SEDs are indistinguishable over 3 - 500 /im; 
the only difference is seen in the near-IR, where the z ~ 1 
SF SED is a factor of a few higher than the z ~ 2 SF 
SED. This region of the SED is probing the stellar bump 
and is a rough proxy for stellar mass. When we compare 
the average sSFR = SFR/Af* of our z ~ 1 and z ~ 2 SF 
galaxies, we find an increase of a factor of ~ 2.5, consis- 
tent with the observed increase in the sSFR with rcdshift 
(e.g., Elbaz et al.||2011 ). In addition, galaxies with more 
dust emission (as evidenced by a higher IR luminosity) 
will have more extinction of the stellar bump. Beyond 
these two effects, there is no evidence from our analysis 
that the SEDs of IR luminous SF galaxies evolve between 
z = 1 and 2, though we are comparing LIRGs at z ~ 1 
with ULIRGs at z ~ 2. The lack of any strong differ- 
ences between our z ~ 1 and z ~ 2 SF SEDs is especially 
interesting, since the former are primarily normal main 
sequence galaxies, while the latter contain a mix of main 
sequence and starburst galaxies according to the IR8 pa- 
rameter. This suggests that mergers either might not 
significantly alter the distribution of dust and gas in high 
redshift ULIRGs, or might no t be as prevalent as they 
are in local ULIRGs. Recently, Diaz-Santos et al. (2011 1 
analyzed the extended mid-IR emission tor a large sam- 
ple of local LIRGs and ULIRGs, and found that once the 
central component of the galaxy is removed, the average 
SEDs of LIRGs and ULIRGs are all rema rkably similar 
in sha pe to the local starburst template of |Brandl et al.| 
(2006). Their results imply that the differences between 
LTRG and ULIRG SEDs are due to processes taking place 
in the cores. High redshift ULIRGs are more extended 
than their local counterparts, so the extended emission 
could be dominating the shape of the SED, causing little 
difference between the z ~ 1 and z ~ 2 SEDs. 

On the other hand, we do see an evolution in the SEDs 
of SF galaxies when we compare our high redshift com- 
posite SEDs to those of local ULIRGs. The general trend 
is that we find more of the IR luminosity is coming from 
colder dust at high redshift, whereas local ULIRGs have 
more shorter wavelength far-IR emission from warmer 
dust. This was originally observe d for high redshift sub- 
millimeter selected galaxies (e.g., Pope et aT||2006[ ), and 
we now show that it extends to a larger high redshift 
ULIRG population selected in the mid-IR. We argue that 
this result would not be as robust without modeling the 
IR emission with multiple dust components since a sin- 
gle dust component would bias the far-IR peak to slightly 
shorter wavelengths, and would fail to properly account 
for all of the substantial warmer dust. The presence of 
more cold dust in high redshift ULIRGs is likely to be in 
part because they are more extended than local ULIRGs 
(even major mergers at high redshift are more extended) 
allowing for a larger fraction of the dust to remain at 
cooler temperatur es ( Daddi et al.|2010 Engel et al.|2010 
Ivison et al. 2011). An interesting test ot our empirica 



SED templates would be to compare them to the ex- 
pected SEDs from radiative transfer models that take 
into account the larger sizes of high redshift ULIRGs. 

6. SUMMARY 

We have combined deep photometry from 3.6 to 
500 /im with Spitzer mid-IR spectroscopy for 151 high 
redshift (z > 0.5) (U)LIRGs in order to explore the rela- 
tive contribution of AGN and SF activity to the total IR 



18 



emission. Our sources bright in the mid-IR (93% have 
S*24 > 100 /xJy), and we have shown that they are rep- 
resentative of the full 24 population above this flux 
limit. We perform spectral decomposition on the mid-IR 
spectrum of each source to determine whether the mid- 
IR luminosity is dominated by continuum emission from 
an AGN or PAH emission from ongoing star formation. 

Based on the mid-IR spectral decomposition, we divide 
our sample into 4 sub-samples (z ~ 1 SF, z ~ 2 SF, sil- 
icate AGN, featureless AGN) and use 95 sources with 
complete spectral coverage to create composite SEDs 
over the full IR range for each sub-sample. We fit a 
two temperature component modified blackbody model 
to our composite SEDs and calculate the dust tempera- 
tures and IR luminosities. By comparing the properties 
of the SEDs, we conclude the following: 

(i) We find no significant difference in the SEDs (PAH 
features, dust temperatures) between SF systems 
with similar stellar masses at z ~ 1 (average Lir = 
4.2 x 10 n L Q ) and z — 2 (average L m = 2.0 x 
1O 12 L ). 

(ii) AGN dominated galaxies with a pure power law 
spectrum are dominated by a warm (~ 100 K) dust 
component in the IR, whereas AGN with silicate ab- 
sorption have a strong warm dust component and 
show a significant cold dust component. 

(iii) We find warmer (~ 20 K) effective average dust 
temperatures for our AGN SEDs than for our SF 
SEDs. 

(iv) We find that a single temperature model is a poor 
fit to the data and neglects significant emission 
at shorter far-IR wavelengths leading to inaccurate 
Lirs and SFRs. 

We compare our composite SEDs created from high 
redshift (U)LIRGs with local te mplates, including the li- 
brary of Chary & Elbaz (2001), and find that the local 



templates differ significantly from our SEDs, particularly 
at far-IR wavelengths where ULIRGs at high redshift 
contain a higher fraction of cool dust than local ULIRGs. 
The difference in warm dust emission between local and 
high redshift SEDs could be attributed to different evo- 
lutionary paths. As the local SEDs used in this work are 
extrapolated from 160 - 850 /im, far-IR Herschel pho- 
tometry is needed for stronger comparisons to be made 
between the cold dust emission in local and high redshift 
galaxies 



(2011) 



We use the criteria presented in Elbaz et al. 
to assess the 'starburstiness' of our sub-samples and 
find, according to either the IR8 parameter or sSFR-z 
relation, that most of our sources lie on the main 
sequence of star formation. The properties of our SF 
and AGN sub-samples are consistent with scenarios 
which link the evolution of star formation and rapid 
black hole growth in galaxies at high redshift. 
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APPENDIX 

A. FAR-IR PHOTOMETRY OF OUR IRS SOURCES 

In the table below, we present mid-IR and far-IR photometry for our 151 IRS sources. We list MIPS 24 and 
70 fxm, PACS 100 and 160 ^m, and SPIRE 250, 350, 500 /xm photometry. We list the results of the mid-IR spectral 
decomposition for each source as a percentage of the mid-IR luminosity that is due to the power-law component. We 
also list the redshifts we calculated from the IRS spectrum. Finally, we state which subsample each source belongs to, 
and we indicate the sources that were not used in the creation of the composite templates. 
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B. CORRELATIONS BETWEEN MODIFIED BLACKBODY PARAMETERS 

In the plots below, we show the correlation between parameters when fitting our two-temperature modified 
blackbody model from Equation (1). We show the results from our Monte Carlo simulations (see s£3j> for the z ~ 2 SF 
galaxy subsample as an illustrative case; the results are similar for the other three subsamples. We overplot the mean 
value of each parameter as the red cross. The plots below show that certain parameters are highly correlated (e.g., 

the temperatures of the warm and cold components) while others are not (e.g., the cold temperature and 
normalization). The errors we report for the dust temperatures and IR luminosities include the covariance between 
the parameters but we note that the error is dominated by the diagonal elements of the covariance matrix. The 
warm and cold temperatures, though correlated, do not overlap in their respective range of values, giving us 
confidence that we are clearly measuring two separate temperature components. 
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